Localized gap soliton trains of Bose-Einstein condensates in an optical lattice 
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We develop a systematic analytical approach to study the linear and nonlinear solitary excitations 
of quasi-one-dimensional Bose-Einstein condensates trapped in an optical lattice. For the linear case, 
the Bloch wave in the nth energy band is a linear superposition of Mathieu's functions ce„_i and se„; 
and the Bloch wave in the nth band gap is a linear superposition of ce„ and se„. For the nonlinear 
case, only solitons inside the band gaps are likely to be generated and there are two types of solitons 
- fundamental solitons (which is a localized and stable state) and sub-fundamental solitons (which 
is a lacalized but unstable state). In addition, we find that the pinning position and the amplitude 
of the fundamental soliton in the lattice can be controlled by adjusting both the lattice depth and 
spacing. Our numerical results on fundamental solitons are in quantitative agreement with those of 
the experimental observation [Phys. Rev. Lett. 92, 230401 (2004)]. Furthermore, we predict that 
a localized gap soliton train consisting of several fundamental solitons can be realized by increasing 
the length of the condensate in currently experimental conditions. 

PACS numbers: 05.45.Yv, 03.75.Kk, 03.65.Db 



I. INTRODUCTION 

Loading Bose-Einstein condensates (BECs) in an opti- 
cal lattice formed by a laser standing wave has received 
increasing interest in the study of nonlinear atomic optics 
Ul) @) S Hi ■ Understanding the properties of BEC in an 
optical lattice is of fundamental importance for develop- 
ing novel application of quantum mechanics such as atom 
lasers and atom interferometers [5|, la, |7|, [sj, |9|, |lO|, lll|. 
Theoretically, some approximation methods are bor- 
rowed from solid state physics, which are used to inves- 
tigate the dynamics of this system. It is mainly due to 
the fact that there are considerable resemblances between 
BEC droplet localized in an optical lattice and electron in 
a lattice. According to the theory of solid state physics, 
there exist band gaps between adjacent energy bands in 
the band structure of solid. In general, the energy bands 
exhibits spatially oscillating phenomena. As discussed in 
Refs. [1211 and [13l|, however, it is possible to generate 
soliton in the band gap when the nonlinearity compen- 
sates for atom dispersion caused by inter-site tunneling. 
The band gap soliton can be called gap soliton. The 
existence of the gap solitons was first predicted based 
on coupled- mode theory [1J|, in analogy to optical gap 
solitons in Bragg gratings. Such a prediction was vali- 
dated by a number of groups using some approximation 
approaches, such as tight binding approximation [ij], a 
complete set of on-site Wannier states [15|, an effective 
mass formula j,16j and plane wave method 113]. Although 
they provide a convenient way to study the gap soliton 
of the BEC, the validity depends greatly on the nature 
of the underlying problem. From this point of view, it 
is desirable to develop a method that does not rely on 
above approximations flSj . 



Strictly speaking, an accurate solution can be obtained 
by exactly solving the full nonlinear Schrodinger equa- 
tion with a periodic potential. However, it is very diffi- 
cult to derive analytical solutions because the full non- 
linear Schrodinger equation is nonintegrable [19| . Conse- 
quently, some asymptotic approaches and numerical sim- 
ulations are used to investigate this question. Using mul- 
tiple scale method, Konotop and Salerno [l8| predicted 
that bright solitons could come into being in a BEC with 
a positive scattering length and dark solitons could be 
stable with a negative scattering length. Subsequently, 
these predictions were proved by using asymptotic theo- 
ries J20|, |21| . Employing numerical simulations, Louis et 
al. [22] analyzed the existence and stability of spatially 
extended (Bloch type) and localized states of a conden- 
sate in the band gaps of the linear Bloch- wave spectrum. 

Especially, Eiermann et al. |2i| reported that gap soli- 
ton do neither move nor change their shape and atom 
numbers during propagation. That is to say, the gap 
soliton is pinned in an optical lattice without attenuation 
and change in shape. Such a soliton can be regarded as a 
spatially localized gap soliton, which is also called funda- 
mental soliton in Ref. 2J]. More importantly, localized 



bright solitons would be very useful for future applica- 
tions, such as atomic interferometry [25] . Subsequently, 
some explanations to this observation were proposed by 
using numerical simulations (see Refs. 
references therein). 
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To better understand the characteristics of the lin- 
ear and nonlinear solitary excitations of quasi-one- 
dimensional (ID) BECs trapped in an optical lattice, we 
develop a multiple scale method to derive analytically an 
explicit expression of the wave function. It is found that 
there are two types of gap solitons in the band gaps. One 



is fundamental soliton, which is always stable and pins 
a fixed position; the other is always unstable and decays 
gradually due to losing a part of its atoms. The paper is 
organized as follows: In Sec. II, we derive ID amplitude 
and phase equations from the original three-dimensional 
Gross-Pitaevskii (GP) equation. Subsequently, by ana- 
lyzing the stability regions of soliton formation, we ob- 
tain the formation condition of the fundamental solitons 
in the band gaps. A novel linear dispersion relation aris- 
ing from the ground state and sound speed of this system 
in the band gaps are obtained in Sec. III. In Sec. IV, we 
develop a multiple scale method to study the nonlinear 
dynamics of the system. We derive a solution of the wave 
function and discuss its dynamical stability in the band 
gaps. It is found that the pinning position and ampli- 
tude of the fundamental solitons are controlled by ad- 
justing both lattice depth and spacing. Furthermore, we 
propose an experimental protocol to observe a localized 
gap train consisting of several fundamental solitons in the 
condensate under currently experimental conditions. A 
brief summary is given in section V. 



II. EQUATIONS OF AMPLITUDE AND PHASE 

Based on mean-field approximation, the time- 
dependent GP equation of full BEG dynamics reads 



where Vq = Ei^/{hjj\^, D — a^^d, fi = Wz/wj, <C 1 and 
Q = 47ras/ao. Expressing the order parameter in terms 
of modulus and phase, i.e., -0 = -ynexp(i$), and then 
separating real and imaginary parts, we obtain 
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Equations @ and ([5]) are (3-|-l)-dimensional, nonlin- 
ear and dispersive equations with a variable coefhcient. 
To solve these equations, we introduce some reasonable 
approximations. Gonsidering a *^Rb condensate in a 
cigar-shaped trap with the frequencies of lJz = 27r x 0.5 
Hz and cj± = 27r x 85 Hz [13], we get n « 0.006. The 
value of Vt is so small that the variation of the pro- 
file of the order parameter is slow in the z direction. 
Thus, the wave function can be separated by -(/'(r, z, t) = 
Ga{r)(l){z,t) with <j){z,t) = A{z,t)eyi.'p[—iijd + iip(z,t)]. 
Here, the modulus and the phase are ^/n = GQ{r)A{z,t) 
and $ — —jit + (p{z,t) respectively. Owing to the 
strong confinement in the transverse direction, the spa- 
tial structure of function Go (r) can be well described by 
a solution of two-dimensional radial symmetric quantum 
t-2 harmonic-oscillator equation, i.e., V^Gq -I- 2Go — r^Ga — 

[--— V^-|-y(X,y, Z)+g|*p]*, (1) 0. The ground-state solution has the form Go (r) = 
2m 



where '^{X, Y, Z, T) is the order parameter of condensate, 
and {X, Y) and Z are the directions of strong transverse 
confinement and lattice. N = J dr\'^\^ is the total num- 
ber of atoms, and g = Airh'^as/'m is inter-atomic inter- 
action strength with the atomic mass m and the s-wave 
scattering length ag {ug > represents the repulsive in- 
teraction). The combined potential V{X,Y, Z) of the 
optical lattice and magnetic trap is 



Gexp(— r /2), where G = I/x/tt can be found from the 
normalization condition J_ jGoprdr — 1. Substituting 
them into Eqs. HJ and ([5]), we obtain 
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where i?^ = X^ -I- F^, Eq is the lattice depth, d = Al/2 
is the lattice spacing, where Al is the wavelength of laser 
beams, ujz and wj, are frequencies of the magnetic trap 
in the axial {Z) and transverse (X and Y) directions, 
respectively. By introducing the dimensionless variables 
t = oj±T, (r, z) = Oq {R, Z) with transverse harmonic 
oscillator length a^ = y^h/mu!±, and ■0 = \/ox>/N^i we 
obtain the following dimensionless GP equation 

.9-0 1„9 , rx. ■ ■?r^^Z^ 



with Q' — Q/{2tt) — 2as/aa. In order to obtain ID 
amplitude and phase equations, we have multiplied Eq. 
dll) by Gq and then integrated the resulting equation once 
with respect to the transverse coordinate to eliminate the 
dependence on transverse plane. An approach similar to 
this one has been widely used in quasi-lD (cigar-shaped) 
BEG problems ^lUliO]. 
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III. LINEAR BLOCK MODES 

We now consider a BEG trapped in a ID optical lat- 
tice. Due to the strong confinement in the transverse 



direction, the system is similar to a wave guide, in which 
the excitation propagates in the elongated direction 31 1 . 
The strong confinement also ensures the dynamical sta- 
bility of the linear excitation 31]. Therefore, we set 
A = uq{z) + a{z, t) with a{z, t) — ao exp(iQ;i) + c.c. and 
(f = ipoexp{iai) + c.c. Here, c.c. is complex conjugate 
and ai = kz — Lot. k is the wave number and w is the 
eigenfrequency. Without loss of generality, we assume 
uo{z) characterizing the condensate background. Con- 
sidering that ao and ipa are small constants, we obtain 
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from the linearization of Eqs. ([6]) and ([T]). Under the 
linear cas e, Q ' « 0. Equation ([9]) is turned into Mathieu's 
equation [32, l33[ 
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with ?7 = TTz/D, q == -VqD'^ /{2'k'^) a,nd p ^ q + 2D'^{fi - 
l)/7r^. Based on the Floquet-Bloch theorem, uq can be 
represented a superposition of Bloch waves, i.e., UQ{rj) = 
&i exp(w77)uoi(?7)-|-&2 exp(— ii^?7)uo2(?7), where woi(?7) and 
■"02 (^) are Mathieu's functions (ce„ or se„), 6i and &2 
are arbitrary constants, and z^ is a Floquet exponent. If 
cos(r/7r) = ±1, the solutions of Mathieu's equation are 
periodic functions and can be expanded as Fourier series 
(detailed expression in Ref. [33'). 

It should be mentioned that in recent experiments the 
characteristic lattice spacing d is determined by the angle 
between the intersecting laser beams forming the lattice 
and varies in the range 0.4— 1.6^to [sJI . The lattice depth 
Eq scales linearly with the light intensity, and varies be- 



tween and E^ 
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is the lattice recoil energy [3J]. So, the dimensionless pa- 
rameters Vq and D are in the range of < Vb < 7.0 x IQ-^ 
and 0.3 < D < 4.05. For convenience, we here set 
D = 3.14 in our calculation. 

From the eigenvalues p and q of Mathieu's function 
(the text-book analysis can be found, e.g., in Ref. [32j). 
figure 1 presents the band-gap diagram for the extended 
solutions of Eq. (|lip which describe noninteracting con- 
densed atoms in an optical lattice. The results are pre- 
sented for the parameter domain (Vo,/i) relevant to our 
problem. One can find that the energy bands (shaded 
areas) are separated by the band gap regions. In these 
band gaps, unbounded solutions exist. The band edges 



FIG. 1: (Color online). Bloch band of BEC in an optical 
lattice in the linear regime is functions of the optical depth Vb 
and the chemical potential fi. The shaded areas denote the en- 
ergy bands. The regions 1 and 2 represent the first and second 
energy band respectively. The areas / and // represent the 
lowest two band gaps in the spectrum. Band edges (solid line) 
a, b, c, d, e, /, ■ ■ • correspond to the eigenvalues and eigenstates 
Po"(ceo), Pr(sei), Pf (cei), P|(se2), P|(ce2), PUses), ■■■ of 
Mathieu's equation, respectively. 



correspond to exactly periodic solutions. The regions 1 
and 2 represent the first and second energy band respec- 
tively, while the regions / and // denote the first and 
second band gap respectively. From Fig. 1, we can con- 
clude that the Bloch wave uq in the nth energy band is 
the linear superposition of the Mathieu's functions ce„_i 
and se„, and uq in the nth band gap is the linear super- 
position of cCn and se„. As is shown below, a complete 
band-gap spectrum of the matter waves in an optical lat- 
tice provides important clues on the existence and the 
stability regions of solitons. 

We next consider the case of Q' ^ 0, and discuss the 
stability problem of soliton formation. Utilizing Eqs. ([5])- 
Pll)l . we obtain 
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Setting Lo — ujr + i^i (where the subscripts de- 
note the real and imaginary parts) [35|, one gets 
that Lol = (7/4)^fc4-h(4fc2/M2)(|a)2 + [e^i/i,) and 

^2 = (7/4)y'A:4 + (4F/m2)(||,)2 _ (fc2^/4), where 7 = 

[l/(2uo)]^ + fcV2 + 2Q'ul. If the imaginary part of 
quasi-particle frequency is a nonzero value, the corre- 
sponding Bloch wave exhibits exponenti al g rowth and 
hence the state tp is dynamical instability 22| • If the fre- 
quency of the associated quasi-particle spectrum is real, 
the soliton would be stable. From the expression of the 



imaginary part of the frequency, one can see the depen- 
dence of the instability growth rate uji on k and J- ^j- . 
On the one hand, when fc = 0, one finds uj ~ 0, which 
is an inessential solution. On the other hand, it is im- 
possible for ^ ^j- being equal to zero because uq is the 

Bloch wave in the energy bands or band gaps. If -^^^ 
is a purely imaginary number (also obtained from Eq. 
(fT2|) ). the wave function ip possesses dynamical stabil- 
ity. We therefore conclude that the stable condition of 
soliton formation is uq = exp(i/?z) where /? is an arbi- 
trary real constant. Because the Bloch wave in the nth 
energy band is the linear superposition of the Mathieu's 
functions ce„_i and se„, it always not satisfies the stable 
condition. Only if Bloch wave in the nth band gap has 
the form of uq — /3ce„ -|- /3isen, soliton possesses dynam- 
ical stability. Thus the linear dispersion relation of the 
nth band gap is 
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Under long-wave approximation, the sound speed is 



Va 



,. du! 

hm — - 

k — ►o dk 



± 



Q'ul 



4 ^d' ' 



(14) 



where the positive (negative) sign represents the right- 
ward (leftward) propagation of the wave packets. For 
the case of D — > oo, the external potential would be 
a harmonic potential [see Eq. ^] and a corresponding 
sound speed is Vgh — \/Q'uq in our notation. This be- 
havior is consistent with both the experimental |36j and 
theoretical results [30|. Obviously, the second term un- 
der the radical sign in Eq. p4|) is arisen from the optical 
lattice potential. The sound speed is the largest in the 
first band gap and gradually decreases with n. Generally 



speaking, the value of 1/D in the experiments [23j, |34| 
would be larger than that of Q' . It implies that the lin- 
ear dispersion relation and sound speed are dependent 
mainly on the lattice spacing. 



IV. NONLINEAR BLOCH MODES 



A. The explicit expression of the wave function 



variables. That is to say, the amplitude and phase 
of the wave function are sought for the forms of A = 
uq{zq,S„t) + e[a(o)(2:o,C,r) -f e'^a^'^\zo,S„T) + • • •] and 

ip = e2[(^(o)(zo,^,r) -t-eV^^H-^cC^T) H ], respectively, 

where the small parameter e represents the relative am- 
plitude of extended states in BEG. Slow variables ^ = 
e{z — Vgt) and r — e^t characterize the slow variation 
of soliton dynamics. Fast variable Z{) — z denotes the 
propagation direction of the lattice wave packets. V^ is a 
group velocity. By substituting them into Eqs. ([6]) and 
([7|), and then separating them in terms of e, Eq. ^ can 
be written as 
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To better understand the nonlinear dynamics of BEG 
in an optical lattice, we here develop a multiple scale 
method to derive an explicit expression of the wave 
function of the condensates in an optical lattice. By 
means of asymptotic expansion in nonlinear perturba- 
tion theory, we propose that the amplitude and phase 
can be expanded by multiple scale methods. In the 
case of that, mathematically, any parameter can be de- 
fined as a function of fast and slow variables, we pro- 
pose each order parameter of the amplitude and phase 
can be written to a function of a fast and two slow 
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From Eq. (|T5|) . one can see that uq is indepen- 
dent on ^. Due to the fast and slow varies possess 
different physical connotation, so it is reasonable that 
these order parameters are written to arithmetic multi- 
ply of function of the fast and slow variables. We may 
set uq{zq,t) = uoi{zq)uo3{t). Similarly, a'-'^' and ip^"^' 
are the forms of a^*''(zo,^, r) = aii{zQ)ai2{S,)ai3{T) and 
(p(')(zo,C,T") = ^ti{zo)Vi2iO'Pi3{T), respectively, where 
i — 0,1,2 ■ ■ ■. Note that the form of Eq. (fT9|) is the same 
as that of Eq. ([9]). In view of the fact that BEC in the ex- 
periments are dilute and weakly interacting: njosl^ ^ 1 
where n is the average density of the condensate, so Eq, 
(fT9| can also be transformed into Mathieu's equation un- 
der the consideration of weak nonlinearity. The solutions 
of Mathieu's equation have been discussed in section III 
By comparing Eq. ^ with Eq. (gOl), one finds a(°) = 



From Eq. ^6 
(fTZ)) becomes 



we obtain uq — ( 



_ (dv'-' 
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So, Eq, 



where the Bloch waves mqi and group velocity Vg can 
be given in Sec. III. fi is the chemical potential. Equation 
([27| is just an explicit expression of the wave function for 
ID BEC trapped in an optical lattice. 

As is discussed in Sec. Ill, the Bloch wave in the nth 
energy band is the linear superposition of the Mathieu's 
functions ce„_i and sen, which is always not satisfy the 
stable condition of soliton formation. Therefore, the con- 
densate in the energy band region can not generate soli- 
ton, only the condensates in the band gaps may be occur 
the soliton. In the following we discuss soliton dynami- 
cal stabilities of the condensates in the band gaps. The 
stability of soliton in nonlinear systems is an important 
issue, since only dynamically stable modes are likely to 
be generated and observed in experiments. 



B. Soliton properties in the band gaps 



dtpoa 
dr 
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The left hand side of Eq. (p4|) is the differentia- 
tion of (/3o3 with respect to r, while the right hand 
side of Eq. ([M|) is differentiation of tpoi {'■Pm) with 
respect to zq {£_). Obviously, both sides of the equa- 
tion must be equal to a constant A, i.e., ipos = 1/(At), 
and <po2 = -(l/2)AuoiC[t^/p^f^]-i. So, ^(«) = 



[uoi^/{2t)][^]-K Correspondingly, we get 

duoaJT) _ uo3{t) 
dr ~ 2t ' 
Similarly, from Eqs. ^, ^ and 
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Under the transformations zq ^ z, ^ — e{z — 

Vgt), and T — e^t, the perturbation parameters 

can be written as uo(z,t) ~ uoi{z)\/t 

-[{z - y^i)uo/(2t)][^]-i, and 
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Vgi)^l\2^t\^ + \Q'ul)\}^\u,,\\ - ^,o^(|f )-^]n, 
where uq\{^z) is the nonlinear Bloch wave of BEC in an 
optical lattice. Finally, the solution of the dimensionless 
GP equation ([3]) is given by 
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To link our analytical results to real experiments, we 
estimate the values of the dimensionless parameters in 
Eq. ([27|) according to actual physical quantities. We 
consider a cigar-shaped ®^Rb condensate (atomic mass 
1.4 X 10^25 i^g a^jjj^^ ^}jg scattering length 5.3 nm) con- 
taining N — 900 atoms in a trap with Wz = 27r x 0.5 
Hz and w j = 27r x 85 Hz (The data are from the ex- 
periment ji^). The parameter fi « 0.006 < 1 in Eq. 
([3|). It implies that the condensate may be regarded as 
a quasi- ID optical lattice in the direction of a weak con- 
finement. Hereafter the radial radius is determined by 
r = 0.01. Based on that ID optical lattice is created from 
one pair of counter-propagating laser beams in real exper- 
iments, the lattice depth and the lattice spacing depend 
on the peak intensity and the angle of the two identical 
counter-propagating laser beams, respectively. For the 
wavelength A_l = 783 nm, used in Ref. [23| , this angle be- 
tween counter-propagating laser beams would be equal to 
2.3°. The periodic potential is V = Eq sin^ (nZ/d) where 
Eq — 0.70 Erec and d ^ Al/2 are the lattice depth and 
periodicity, respectively. Accordingly, the time and space 
units correspond to 0.2 ms and 3 fim, respectively. These 
units remain valid for other values of TV, as one may vary 
Vq accordingly; in this case, other quantities, such as D, 
also change. Based on these proposed, we obtain the di- 
mensionless parameters D « 0.33, Q' « 9 x lO^'' and 
Vo ~ 7.9. 

As is discussed above, the nonlinear Bloch waves in 
the nth band gap can be given by the linear superposi- 
tion of ce„ and se„ under the case of weak nonlinearity. 
On the basis of the fact that the coefficients of periodic 
Mathieu's functions depend on eigenvalue q [33| (i.e., Vo 
in our notation), we presuppose that the Mahtieu's func- 

1000 1000 

tions arc cei = J2 ^o cos{kri) and sei ~ J2 ^o ^^{kri) 

fe=0 k=0 

in the following calculation. 

First, we choose a linear superposition form of mqi (^) = 





FIG. 2: (Color online). The space-time evolution of the den- 
sity of the condensates in the first band gap. The Bloch wave 
is chosen as uoi{z) = cei + sei. The parameters used are 
the lattice depth Vo = 7.9, interatomic interaction strength 
Q' = 9 X 10"'^, the lattice spacing D — 0.33, the radial radius 
r = 0.01, and the chemical potential \i, — 1.12. All parameters 
are in dimensionless units. 



cei -|- sei as the nonlinear Bloch wave in the first band 
gap. The time-evolution of the density distribution of the 
condensates in this case is plotted in Fig. 2. We see that 
the peak of wave packets decreases exponentially with 
time and eventually vanishes. Thus the state is always 
unstable and be called sub-fundamental soliton in Ref. 
2J| . Such a soliton with a very small initial total number 
of atoms loses a part of atoms with time going on, so it 
is unstable (refer to Ref. |24|1. 

Secondly, we choose uoi{z) = cei +«sei, which satisfies 
the stable condition of soliton formation, as the nonlin- 
ear Bloch wave in the first band gap. Figure 3 shows 
the space-time evolution of the density the condensates 
in first band gap. A strong peak appears in the conden- 
sate with a dimensionless lengths of 12.6 (about 38 lattice 
sites), and maintains its shape and magnitude. It implies 
the existence of a bright gap soliton. As the time going 
on, the bright gap soliton is pinned in the optical lat- 
tice without both attenuation and change in shape. This 
behavior indicates that it is a spatially localized bright 
gap soliton, which is arisen from the interplay between 
the tunneling of periodical potential and nonlinear inter- 
action of the system. Moreover, the width of the peak 
in the z — t plane is found to be a dimensionless length 
with 2.0, i.e., ~ 6 fim in real space. The value is in 
good agreement with that of the experiment observed by 
Eiermann et al. [23|. The agreement illustrates that our 
method can describe the dynamics of BEC trapped in 
an optical lattice very well. The type of bright gap soli- 



FIG. 3: (Color online). The space-time evolution of the den- 
sity of the condensate with 12.6 dimensionless lengths (about 
38 lattice sites), which corresponds to about 38pm in real 
space. The Bloch wave is chosen as uoi{z) — cei -I- isei. 
Other parameters used are the same as Fig. 2 



tons are called fundamental solitons in Ref. [2j] . Similar 
phenomena can also be obtained inside the other band 
gap. 

From the results discussed above, the solitons residing 
in the band gaps are the fundamental or sub- fundamental 
solitons depending on their position in the band gap 
whether the stable condition of soliton formation can be 
satisfied or not. 

In real experiment, the ID optical lattice is created 
from one pair of counter-propagating laser beams, and 
the lattice depth depends on the peak intensity of the 
two identical counter-propagating laser beams. That is 
to say, the lattice depth can be adjusted by varying the 
intensity of the counter-propagating laser beams. We 
here depict how the lattice depth infiuences the funda- 
mental soliton in Fig. 4 (with all the rest of the system 
parameters fixed). From solid line (Vb = 7.9) and dashed 
line (Vq = 10), one sees the amplitude of the fundamen- 
tal soliton increasing with the increasing of the lattice 
depth. Due to both wave packets containing the same 
total number of atoms but the bosons in deeper well are 
captured more tightly, the tunneling probability varies 
smaller [l9|, ISTj. To balance the same nonlinear effect 
of the system, the tunneling rate of bosons in deeper 
well becomes much larger to achieve the same dispersion 
effect, which results in the amplitude of the fundamen- 
tal soliton increasing. Therefore, the amplitude of the 
localized gap soliton increases with the increasing the in- 
tensity of the counter-propagating lasers beams in the 
experiments. 

Subsequently, we observe the soliton characteristics in 




FIG. 4: (Color online). The distribution of the density of the 
condensates with different lattice depth at t = I, where the 
dimensional parameter Vq = 7.9 corresponds to the lattice 
depth in the experiment [23|. Other parameters used are the 
same as Fig. 2. 



a longer condensate. Our numerical calculations are per- 
formed for the condensate cloud in the ground state ex- 
tending over 35 dimensionless lengths (about 106 lat- 
tice sites), which corresponds to 105 /zm in real space. 
The condensate cloud contains about 2.5 x 10^ atoms 
under the consideration that the atomic density keeps 
unchanged. Fig. 5 shows the space-time evolution of the 
density of the condensates in this case. It is shown that 
there exhibits a localized gap soliton train consisting of 
several fundamental solitons in the condensate. Similarly 
to the property of a single fundamental soliton, the soli- 
tonlike wave packets in the train are immobile (i.e., have 
zero group velocity). In reality, the condensate is loaded 
into an optical lattice from a crossed optical dipole trap 
[23|, which results in an initial state of a finite extent. 
A BEC wave packet centered around a particular quasi- 
momentum in a given band is created by ramping up of 
a static lattice with subsequent linear acceleration to a 
given velocity [23|. In order to observe a localized gap 
soliton train, we may propose experimental protocols ac- 
cording to the ex per imental observation of a single fun- 
damental soliton [23|. First, the atoms are initially pre- 
cooled in a magnetic time-orbiting potential trap using 
the standard technique of forced evaporation leading to a 
phase space density of ~ 0.03. Subsequently, the atomic 
ensemble is adiabatically transferred into a crossed light 
beam dipole trap, where further forced evaporation is 
achieved by lowering the light intensity in the trapping 
light beams. With this approach, one can generate pure 
condensates with typically 9 x 10"* atoms. By further 



lowering the light intensity, one can reliably produce co- 
herent wave packets of 9000 atoms. For this atom num- 
ber no gap solitons have been observed. Therefore, one 
removes atoms by Bragg scattering. This method splits 
the condensate coherently leaving an initial wave packet 
with 2.5 X 10'^ atoms at rest. Then, the wave packet 
centered on a particular position in a given band gap 
(which satisfies the stable condition of soliton formation) 
is created by switching off one dipole trap beam from 
a crossed optical dipole trap, releasing the atomic en- 
semble into ID horizontal waveguide with transverse and 
longitudinal trapping frequencies uj±_ = 27r x 85 Hz, and 
t^ll = 27r X 0.5 Hz, and then accelerating the periodic po- 
tential to the recoil velocity Vr = h/mX. This is done by 
introducing an increasing frequency difference between 
the two laser beams, creating the optical lattice. The ac- 
celeration is adiabatic, which results in an initial state of 
an about 105 /xm extent. In view of the fact that the tun- 
neling rate of about 900 atoms extending the length of 38 
fim (in our above simulation) can balance its nonlinear 
energy, such a system generates a fundamental soliton. 
With both the length of condensate and the total num- 
ber of atoms increase, the wave packet exhibits violent 
dynamics. During this evolution the wave packet con- 
taining 2.5 X lO'^ atoms is separated from the surround- 
ing atomic cloud into several BEC wave packets, so the 
periodic structure of a train of the localized wave packets 
emerges. Such a structure represents a train consisting of 
several fundamental solitons, which is supported by the 
combined action of the repulsive nonlinearity and anoma- 
lous diffraction caused by intersite tunneling in the band 
gaps[il,[ii,|33. 

Finally, we study how the lattice spacing influences the 
fundamental soliton or soliton trains as shown in Fig. 
6. In practice the variation of lattice spacing is easy 
to control by adjusting the angle between two counter- 
propagating laser beams. From solid line {D = 0.33) 
and dashed line {D = 0.99), we find that when lattice 
spacing occurs a slight difference, the pinning position 
and the amplitude of each fundamental soliton have a 
little change keeping the distance between adjacent soli- 
tons unvaried. It illuminates that the condensate cloud is 
separated into three wave packets for both D — 0.33 and 
D — 0.99. When the lattice spacing D varies from 0.33 
to 0.99, the condensate from 91 decreases to 30 wells. 
Owing to the center position of each wave packet float- 
ing, the pinning position of the fundamental soliton are 
set to move. And each BEC wave packet still forms a 
localized gap soliton, which comes from the balance be- 
tween the nonlinearity and atom dispersion caused by 
intersite tunneling. However, when the lattice spacing D 
varies from 0.33 to 0.99, the number of atoms in a well 
increases from 28 to 83. For the same length conden- 
sate, the tunneling between adjacent well varies easier 
with the increasing of the atomic number confined in a 
well. To balance the same nonlinear energy of the sys- 
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FIG. 5: (Color online). The space-time evolution of the den- 
sity of the condensate with 35 dimensionless lengths (about 
106 lattice sites), which corresponds to about 105 fim in real 
space. Other parameters used are the same as Fig. 2. 



FIG. 6: (Color online). The distribution of the density of 
the condensates with different lattice spacing D at t = 1, 
where the dimensional parameter D = 0.33 corresponds to 
the lattice spacing in the experiment [23]. Other parameters 
used are the same as Fig. 2. 



tern, bosons gathering around several wells vary easier, 
which results in the amplitude of each localized soliton 
having a increasing trend. When lattice spacing varies 
much larger [see dotted line D — 16.5 in Fig. 6], one 
find that there appears only a fundamental soliton in the 
entire condensate with about 30 dimensionless lengths. 
This mainly reason is that the condensate in this case 
only contains about 2 wells. The tunneling of bosons in 
adjacent lattice achieves the dispersion effect to balance 
the nonlinear energy of the system, so there exhibits only 
a fundamental soliton in this case. 

From the results discussed above, we can conclude that 
the condensate generating a single fundamental soliton or 
a localized gap soliton train consisting of several funda- 
mental solitons can be controlled by adjusting the length 
of condensate or (and) the lattice spacing. Our theoret- 
ical results reported here is important in understanding 
the fundamental soliton physics of BEG in the future. 



V. CONCLUSION 

In summary, we develop the multiple scale method to 
study the linear and nonlinear solitary excitations for ID 
BEG confined in an optical lattice. After averaging over 
the transverse variable, a hydrodynamical model of the 
amplitude and phase is derived. In the linear case, the 
Bloch wave in the nth energy band is the linear superpo- 
sition of the Mathieu's functions ce„_i and se„, and the 
Bloch wave in the nth band gap is the linear superposi- 
tion of ce„ and se„. In addition, we find that the stable 



condition of soliton formation is that the Bloch wave uq 
in nth band gap satisfies uq — [3cen + i(3sen- Under this 
stable condition, a novel linear dispersion relation and 
sound speed are derived. It is found that the linear dis- 
persion relation and sound speed depend mainly on the 
lattice spacing. 

For the nonlinear case, we derive a solution of the wave 
function of the condensates with weakly inter-atomic in- 
teraction, and discuss its stability for condensate ^^ Kb 
in band gaps. It shows that there are two types of gap 
solitons in the band gaps. One is fundamental soliton, 
which is always stable and pins a fixed position; the other 
is subfundamental soliton which is always unstable and 
decays gradually due to losing a part of its atoms. Only 
when the Bloch wave in the band gaps satisfies the stable 
condition, the condensates exhibit the fundamental soli- 
tons, otherwise there appears the sub-fundamental soli- 
tons. Furthermore, the pinning position and the ampli- 
tude of the fundamental solitons in the lattice can be 
controlled by varying the lattice depth and spacing. We 
also propose an experimental protocol to observe a local- 
ized gap soliton train consisting of several fundamental 
solitons for BEG trapped in an optical lattice in future 
experiment. 
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